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Abstract 

We consider a thin soft adhesive interphase between dissimilar elastic media. The material 
of the intermediate layer is modelled by elasto-plastic pressure-sensitive constitutive law. An 
asymptotic procedure, together with a novel formulation of the deformation theory of plasticity 
for pressure-sensitive materials, is used in order to derive nonlinear transmission conditions 
for the corresponding imperfect zero-thickness interface. A FEM analysis of the original three- 
phase structure is performed to validate the transmission conditions for the simplified bimaterial 
structure. 

Keywords: Adhesive joint; Pressure-sensitive deformation theory; Imperfect interface; Nonlinear 
transmission conditions 


1 Introduction 

Adhesive joints are widely used in structural engineering applications involving composite materials, 
such as marine jH, aerospace j3| and automotive structures. In these structures, polymeric 
adhesives are used to join different materials such as plastics, metals, fibre reinforced composites 
and others. 

Early models of adhesive joints tended to assume that the adhesive behaves as a linear-elastic 
solid BS- However, many adhesives (e.g. rubber modified epoxies) exhibit large plastic strains 
before failure Q. A number of yield criteria have been proposed to model the plastic behaviour 
of polymers. Among these, the Tresca and von Mises criteria, originally developed to describe the 
yield behaviour of metals, have been also used. However, these criteria cannot accurately predict the 
behaviour of polymeric adhesive specimens under multiaxial loading, as yielding in these materials is 
sensitive to hydrostatic as well as the deviatoric stress 0,3. As a consequence, pressure-dependent 
plasticity theory is more appropriate in this case. In this context, the Drucker-Prager yield criterion 
is widely used for polymeric materials. 
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The numerical FE modelling of thin adhesive interphases is complicated by the very large aspect 
ratio of the intermediate layer, which requires extremely refined meshes. Alternatively, the so-called 
imperfect interface approach can be adopted. It consists of replacing the actual thin interphase 
between the two adherents by a zero thickness imperfect interface. The interface is supplied with 
special transmission conditions, which incorporate information about geometrical and mechanical 
properties of the original thin interphase. 

Although imperfect transmission conditions for elastic interphases have been intensively inves¬ 
tigated BQ, analytical results for elastoplastic adhesive interphase are very scarce and mainly 


limited to rigid adherents and pressure-independent plasticity |13l . Il4 1. 


In this paper, nonlinear transmission conditions for a soft elasto-plastic pressure-dependent in¬ 
terphase are derived by asymptotic techniques. The interphase material is described by a novel 
formulation of the deformation theory for pressure-dependent media. Early models of pressure- 
dependent deformation theory can be found in 15 1 and (l^. However, the model proposed by 151 


only accounts for elastic and plastic branches separately, i.e. the yield condition is not incorporated 
into the generalized constitutive parameters. On the other hand, the model proposed in |l6l | applies 
only to a few specific loading paths. To overcome these limitations, the classical Hencky deforma¬ 
tion theory is generalized in Section [3] to the case of pressure-sensitive materials by including the 
dependence of inelastic deformation on the first stress invariant, in the spirit suggested by Chen [l7| 


and Lubarda |l8l |. The general theory is then specialized to the the case of Drucker-Prager model 
with associated flow rule and isotropic linear hardening. 

The transmission conditions are validated by numerical examples based on accurate finite ele¬ 
ment simulations. The accuracy and efficiency of the proposed approach prove to be very high for 
different monotonic loading conditions. 


2 Problem formulation and transmission conditions for elasto-plastic 
pressure-dependent interphases 

Consider a structure composed by two dissimilar elastic materials joined together by a thin adhesive 
interphase, see Fig. [T] The thickness of the interphase is small in comparison with the characteristic 
size of the body: h = e/i* L, ^ L. Here e is a small positive parameter, e <C 1. 

The adhesive material is assumed to be soft in comparison with the two adherents and may 
exhibit a very general nonlinear constitutive behaviour, including compressible plastic deformations, 
the only assumption being that the material is isotropic. 

Adopting the deformation theory of plasticity, the constitutive laws of the interphase in the 
plastic regime will be described in terms of nonlinear elasticity as 

— ^^kk^ij T , (1) 

where the generalized Lame parameters are functions depending on the deformation within the 
interphase 

A = A(JbJ|), = (2) 

Note that, due to the isotropy assumption, the dependence is written in terms of the two invariants 

~ ^kki J2 ~ (^) 


2 













L 

Figure 1: Bimaterial structure with thin soft adhesive joint. 


where Cij denotes the deviatoric part of the strain tensor. For the sake of keeping the analysis 
simple, the dependence on the third invariant is omitted. However, if necessary, the model can be 
generalized to include the effect of all invariants. 

In the elastic regime, the elastic constants are assumed to be constant 

A = eAo, /i = e/ro, (4) 

where the elastic parameters Aq, fJ-o are comparable with those of the adherent materials, X±, fj,±. 

We assume in our analysis that, for any admissible deformation of the material within the 
interphase, there exists a constant parameter such that 


0 < J|) <u^< 1/2, 

where Jf, Jl) generalized Poisson’s ratio 

A(Jf,J|) 




2[A(Jf, Jl) 4“ j(Ji) J|)] 


(5) 


( 6 ) 


Under such assumption, the following two imperfect transmission conditions linking the values from 
different parts of the interphase are asymptotically satisfied, see [l6| 


|t2l(xi,X3) = 0, 


(7) 


and ^ 

t2{xi,X3) = —A(xi,X3)M(xi,X3), 


( 8 ) 


where t 2 = [<721, 0 ' 22 ,< 723 ]^ is the traction on the interphase surfaces, u = [ui,U 2 ,U 3 ]'^ is the dis¬ 
placement vector, and the jump of a function / across the interphase is denoted by |/]] = f{x 2 = 
h) — f{x 2 = —h), and the tensor A is given by 


/fi 0 0 \ 

A = 0 A 2/1 0 . (9) 

Vo 0 flj 
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The transmission conditions ([71) and (| 8 |) are valid under the following ellipticity conditions 


0 < /i < e/r*, 0 < A + 2/i < eA*, 


( 10 ) 


where /r* and A* are constants comparable in values with those of the adherent materials. In fact, 
we can take = /tq and A* = Aq + 2fiQ. 

The asymptotic procedure allows us also to estimate the components of the strain tensor, which 
to the leading order 0{e~^) read 


Eij = 0(1), 




A J 1 ^; ^22 


du2 

dxo 


( 11 ) 


By application of the generalized Hooke’s law ([T]), we obtain the following relationships 


£i2 — 


(^i2 


b J I5 ^7 ^22 


<722 


2/1' “ x + 2il' 

Finally, the leading terms of the strain tensor invariants can be calculated as 

<722 


— £22 — 


A+ 2/2 


( 12 ) 


(13) 


^2 — 2^22 + ^12 + ^23 — 


a. 


22 


3 (A + 2 / 2)2 4 


+ - -o + 4 -2 • 


(14) 


Note that the material functions A and / 2 , at this point, are not specified and should be defined 
depending on the particular elasto-plastic model adopted for the interphase material. However, all 
the stress components of the traction vector t 2 as well as the components of the strain tensor Sij 
in their leading terms are functions only of the space variables xi and X 3 . This immediately implies 
that, in the case under consideration, the strain invariants are also functions of those two variables. 


Jf = Jf(xi,X3), J| = J|(X1,X3). 


(15) 


The transmission conditions (| 8 |) are then used to write the strain invariants dnD and (1141) in terms 
of components of displacement jump 




76 — 
Jo — 


12/l2 


Kf + TTTT Nf + TTIT Nsl 


16/l2 


16/l2 


(16) 

(17) 


This allows us to formulate the nonlinear transmission conditions for a thin soft nonlinear interphase 
in their final form 


P 2 I — 0, *2 


—a(M)M, 


(18) 


where the matrix-function A is explicitly defined by A([[u]) = A(A( J^, J|),/2( J^, J|)) and the 
relationships (|16p and m- 

Let us recall that the assumptions (1101) should be satisfied to ensure the applicability of the 
transmission conditions. However, even if the estimates (IlOp cannot be checked a priori, the nonlin¬ 
ear transmission conditions (|18l) can be successfully used in the modelling. In this case, a posteriori 
verification of ( 1101 ) will be necessary in order to confirm the validity of the final results. 
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Summarizing, eqs. (1181) constitute the transmission conditions describing a soft imperfect inter¬ 
face in a bimaterial structure. These conditions allows one to replace the thin plastic interphase 
with an imperfect nonlinear interface (of zero thickness) and consider a simplified outer problem 
for the two bonded materials. The thin interphase influences the solution only via the imperfect 
transmission conditions which incorporates the most important information on the properties of the 
adhesive material and the nonlinear deformations developing within the interphase. 

However, in order to use the nonlinear imperfect transmission conditions (1181) one needs to define 
the generalized Lame parameters (|2]) as functions of the strain invariants depending on the mechan¬ 
ical properties of the interphase material. If the parameters are constant, the conditions reduces 
to the thin soft elastic interphase case. In case of elasto-plastic material obeying the von Mises 
criterion, the analysis was presented in |l4| . In the next section, we derive the proper deformation 
theory for a pressure-dependent plastic material, with particular reference to the Drucker-Prager 
criterion. 


3 Pressure-dependent deformation theory 

3.1 Main Assumptions 

In this section, the classical J 2 deformation theory is generalized to the case of pressure-sensitive 
materials, so that the first invariant Ji will also be included in the formulation. We make the 
following assumptions: 

• The material is initially isotropic. 

• The principal axes of the plastic strain tensor are coincident with those of the stress tensor 

aij. 

• The plastic deviatoric strain tensor is proportional to the deviatoric stress tensor Sij. 

In the following, it is convenient to introduce the volumetric and deviatoric decomposition of the 
stress and strain tensor as: 

— ^m^ij T Sij , 

— ^m^ij T ^ij ) 

and the additive composition of the elastic and plastic strains 

= ^ij + + ^m) ^ij + ( 21 ) 

in order to derive the closed-form relationship. According to assumptions (2) and (3), the plastic 
strain can be written as 

^ij — T (22) 

where tpi and 4>2 are functions which represent the hardening behaviour of the material. In the case 
of plastic loading, it holds that (pi ^ and (p 2 > 0, while in the case of elastic loading (pi = (p 2 = 0 
prevails. Note that Eq. (I22p can be derived following Chen’s suggestions fl^ . 


( 19 ) 

1 
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The constitutive relation between stress and elastic strain decomposed in its volumetric and 
deviatoric parts is given by 


2 G’ 




m' 


(23) 


The material is assumed to be elastic compressible, so that 0 < < 1/2 and K < oo. 

Under the given assumptions, the nonlinear stress-strain law for a pressure-dependent plastic 
material becomes 

where the generalized Lame’s coefficients have been introduced: 


4>2) — 


312 + (02 - (t)i)E 


3{1 + 12 + 02-E)(1 — 2i 2 + (piE) 


E, p-ih) = 


E 


2{i +12 + 4>2Ey 


(25) 


It should be noted here that these coefficients coincide in the pure elastic regime (0i = 02 = 0) 
with the elastic Lame’s coefficients. Further relationships for generalized elastic constants can be 
derived in the same manner: 


^( 01 , 02 ) = 


3E 


3 -|- (202 -|- (j)i)E 


3z^ + (02-0i).F - 


E 


3{l - 2 i 2 + (l)iE)' 


(26) 


To define the functions 0i and 02, let us consider the relationships for the first and second 
invariants of the strain and stress tensor, i.e. 


E — f'p 

where the proportionality relation 




can be used to obtain: 

1 


(27) 


— 301 (Tm — 01~ ~ 02'^2- (^3) 

Thus, the scalar functions 0i and 02 can be obtained in the case of a multiaxial stress state as 

je'P / 

.#1 = - 0 , 4^ = Ew- 

In order that these relations can be used in the analysis of the thin plastic layer, the functions 0i 
and 02 must be defined in terms of total strain invariants. From the nonlinear elasticity relationship 
(j24p . it is noted that the first invariants of the stress tensor and the total strain tensor are linearly 
dependent, 

= (3A + 2/i) Jf = 3k Jf, (30) 


whereas the second invariants of the deviatoric tensors are related by 

^ = 2/i /f. 


(31) 
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This, in turns, allows us to evaluate the generalized material parameters as functions of the invariants 


2fl = 




3K=^. 

J! 


Substituting these relations into Eqs. f|25p and 

Jf 1 - 2u 


J? 


E 


4>2 = 


one can define after some algebra: 
1 + 1/ 




E 


(32) 


(33) 


These relationships are more useful in the analysis than ([29]). This also eives the well known 
relationships between the invariants of the plastic strain and the total straiiTJ: 




E 


(34) 


In the next section, we develop the deformation theory in the case of uniaxial stress state. In Section 
13.31 the deformation theory is generalized to the case of general multiaxial state and Drucker-Prager 
material. 


3.2 Uniaxial Stress Test 

The uniaxial stress test is one of the most common tests, widely accessible, and it may provide at 
least two sets of experimental data: ax{Ex) and £y{ex)- The specific form of the scalar functions (j)i 
and 4>2 for the case of a uniaxial stress state can be derived in the following way. Assuming that 
the direction of loading coincides with the x-direction {ux 7 ^ 0; = e§), the scalar functions (()i and 

02 become 

V2eI ^ \el-el\ 


01 = 


02 = 


(35) 


Let us assume in the following an elastic-plastic material with linear hardening as shown in Fig. |2] 
The flow stress ax can be expressed in terms of the plastic modulus E^ as (cf. Fig. [SJ)) 


axiel:) = EP£P + as, 

whereas the transversal deformation can be expressed as (cf. Fig. |2]i) 


(36) 


SyieP) = -uPeP - 


E 


a. 


The plastic part of Sy can be derived as 


£Py = £y-£l = - —as + i 2 —=-\iy^- ly— I e'x, 


re _ _,,Pc-P 


(Tx 


EP 


(37) 


(38) 


^ The last formula can be also justified by the fact that the tensors e®, and e are proportional to each other: 

1 


= 2fx(j>2e'' = 


1 + 2/r02 


e. 
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Figure 2: Example of constitutive behaviour in a uniaxial stress test for a linear hardening material: 
(a) Axial stress ax as a function of the total axial strain Ex', (b) Axial stress ax as a function of 
the plastic axial strain Ex', (c) Transverse contraction Ey as a function of the total axial strain Ex', 
(d) Transverse contraction Ey as a function of the plastic axial strain Ex- Three cases of transverse 
contraction are considered: plastically incompressible material = 0 = 0.48 and = 0.53), 

null plastic transverse contraction Ey = 0 = 0.027 and = 0.03) and an intermediate case 

= 0.2 and = 0.22). 


so that the functions (j)i and (j )2 are obtained in terms of the plastic strain in a uniaxial stress state 
as 


M^x) 


{E - 2vPE + 2 vEP)e^x 
E[EPEx + < 7 ^) 


02 (e?) 


{E + uPE - i'EP)e^x 
E{EPEx + ( 7 s ) 


(39) 


Using the additive composition of the elastic and plastic strains, i. e. Ex =£% + £x, the stress ax 
and the transversal deformation Ex are obtained as functions of the total strain Ex 


where 


E - E^P 

^x{£x) — E ^£x H (7s 

h/ 


£y{£x) = - 


- IJ^P 


-(7s 


E^P 


EEP 
E + EP' 


uPE 
E + EP' 


(40) 

(41) 
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Finally, the dependence of 4>i and 02 upon the total strain is given by 


01 (ex) 

02 (ex) 


l-2vP + 2u- 


EP\ ESr — Co 


E I EPEx + 


-{l + uP-u— 


EP\ Esx — Co 


E I EP£x + Cs 


(42) 

(43) 


The functions 0i(e3,) and 02 (e^,) corresponding to the material parameters presented in Fig. |2]are 
shown in Fig. (3] Interestingly, in the case of uniaxial stress state the ratio of the functions 0i(e2;) 



Figure 3: (a) Function 0i(ex) and (b) function 02 (e^) corresponding to the material parameters 
presented in Fig. [2] 


and 02 (^x) is a constant value: 


01 (ex) _ {l-2vP)E+ 2vEP 
4>2 {£x) {l + vP)E -vEp 


(44) 


3.3 Pressure-dependent yielding: Drucker-Prager Material 


Consider a material obeying the Drucker-Prager yield criterion, 

F( Jf, J|) = aJf + yj| - (45) 


where a > 0 is called pressure-sensitivity index and assumed to be a constant while kg is a positive 
function of the plastic strain 

kg = kg{£^) > 0. (46) 

In the case of uniaxial stress test 

-I- 

Assuming the normality rule (plastic strain increment proportional to the normal of the yield surface 
riij = 2a^fj^5ij + Sij) and proportional loading we have (cf. (l8| l 

eC = 4>2{2a^7^5ij + Sij), (48) 


^s(^x) — 


(E^£^ + Cs) 


(47) 
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so that 


(49) 


Jf — Qa.(p2\/~J^-, y ^2 ~ 

Note that the ratio between the two plastic strain invariants is constant 


Jf 


= 6a. 


(50) 


The general setting presented in Section 13.11 can be specialized to this type of material by defining 


6 ayjf (/>2 

91 = - 7 ^ -) 92 = 92 - 

It is important to emphasize here that both functions (^i and (/>2 equal to zero simultaneously in the 
elastic regime or are both different from zero in the plastic zone. 

We need to represent now the functions (pi and (j )2 depending only on the invariants of the total 
strain tensor. Let us assume that we have been successful with the second function, 



(f>2 = 4>2{JI,J2)^ 


(52) 


then the remaining function (pi can be computed from Eq. (|51l) under taking into account Eqs. 
and dMI 

, f, ^/^p2 2jl 1 — 2z^ + piE 

9i = 6a ---= 6a92 —=- 77— = 6a92' 


or 




Pi=Pi{JlJ-) = 


2>K Ji 1 + zz + P 2 E Ji 

6 a(l - 2v)p2\fJi 


(53) 


(54) 


Jf(1 + v) + (P 2 E{JI - 6ay^)■ 

Using Eqs. (1521) and (|54p . the generalized elastic constants defined in Eqs. (|25l) and (|26p can be 


(55) 

(56) 

(57) 

(58) 


rewritten in terms of the only two invariants of the total strain tensor: 

Jl{p2) = G{p2) = 

E{P2) = 
y{p2) = 


2{l + v + p 2 Ey 
3E Jl{l + V + P 2 E) - 6ap2E^/Jl 
1 + + P 2 E J^(3 + 2 P 2 E) — 12ap2Ei/iy^ ’ 
Jf (3z^ + P 2 E) - 6ap2E^/Jl 
Jf(3 + 2 P 2 E) - 12ap2Ei/J^’ 


~ _ Jf (1 + U + P 2 E) - 6ap2Ey^ 

3Jf(l-2z/)(l + z/ + (/. 2 .E) 


Finally, the multiaxial hardening law for a linear hardening material (cf. Fig. [2|) can be obtained 
from the uniaxial one, Eq. (|47l) . by substituting Ex with the equivalent plastic strain, which can be 
defined for the Drucker-Prager material as 


iP = 


V3E 


(1 + uP)E - vEP 


I j-gp _ 2(\/3 + 3a) 
'h — 


(59) 
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This gives 


( 60 ) 



where 

VSEPE _ 2EP{V3 + 3a) 
~ (1 + vP)E - vEP ~ 3 


(61) 


In case of the uniaxial stress test discussed in the previous section, Eqs. fl44p and m allow us to 
determine the unknown parameter a > 0 from the experiment 


Jf (1 - 2uP)E + 2vEP _ 1 (1 - 2vP)E + 2 vEP 

6 yj| (1 + yP)E - vEP ~ ^ (1 + uP)E - uEP ^ 

Thus, to complete the analysis, the only unknown function ^>2 should be defined. In general, this 
function depends on both strain invariants and J| 


</>2 = </>2(Jf,J|). 


(63) 


To this end, note that we can eliminate Jf and jf from Eqs. fl34p i . (1501) and then solve Eqs. fjM)) i . 
dsni) to obtain J 2 ) '^2^('^i! ’^ 2 )- function (/>2(Jf, J|) is then derived as 





3a(l + i')EJl + (1 — 2z/) \2>EyJlJ^ — (3a + \/3) (1 + iy)as] 

E ^ISa'^EyO^ — 3a[iilJf — (1 — 2v){u-,JlJ^ + fT^)] + \/3(l — 2z^)(a;yO| + cjs)} 

(64) 


Finally, the function i^i(Jf, J|) is obtained from (1541) 




6 a {3a(l + v)EJl + (1 — 2v) — (3a + -\/3) (1 + i2)as] } 

E [3E{Jf — Qa^fj^) + (3a + '/s) (1 + + 6acJs)] 


(65) 


The evaluation of the functions (pi and (j )2 of the two strain invariants Jf and J| for a Drucker- 
Prager material having the uniaxial stress-strain curve given in Fig. ^{E — 813 MPa; u — 0.3; as 
— 50 MPa; EP — 81.3 MPa; u^P = 0.2; uP = 0.22) is shown in Fig. ID It is noted that the function 
(pi shows a discontinuity line emanating from a pure shear state on the yield surface and separating 
regions of different sign. The function (p 2 also shows a discontinuity line emanating from the vertex 
of the Drucker-Prager cone. 

The deformation theory developed in this section for Drucker-Prager material with linear hard¬ 
ening has been validated by comparison with the flow theory for several monotonic loadings. Finite 
element simulations on a single element with displacement control have been performed in ABAQUS 
and results (not reported here) are always in agreement with the deformation theory. 

We present here two examples of monotonic loadings, designed in such a way to investigate 
the discontinuities shown by the functions <pi and (p 2 - It will be seen that these discontinuities are 
genuine features of the Drucker-Prager model with associated flow rule and are obtained also by 
the more rigorous flow theory of plasticity. 

In the first example, the following monotonic loading has been applied in two steps (see Fig. [5]): 
first a dilatational deformation (from point O to point A in Fig. [5^), followed by an isochoric 
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Figure 4: (a) Function <j)i and (b) function <p 2 of Drucker-Prager’s material with linear hardening as 
functions of the two invariants of strain, Jf and J^- The dashed red line indicates the elastic limit. 


deformation (shearing in the plane X 1 -X 2 , from point A to point B in Fig. |5^). This specific loading 
path has been chosen as it crosses the discontinuity line for the function (pi. The value of the 
function <pi along the isochoric path AB is shown in Fig. [S]^: results of finite element simulation 
and deformation theory are in agreement. 




Figure 5: (a) Function cpi in the plane: elastic domain = 0 (light grey), plastic domain 

with a discontinuity line separating the regions <pi > 0 (light blue) and (^1 < 0 (light yellow). The 
red dashed line indicates the loading path, (b) Function (pi along the isochoric path AB. The 
discontinuity is clearly visible and results of finite element simulation and deformation theory are 
in agreement. 


In Fig. [6l we present the stress-strain curves along the isochoric path AB of Fig. |5^. In particular. 
Fig. |6^ shows the shear stress ai 2 and Fig. [G]^ the normal stress an, both as a function of the shear 
strain £12. The agreement between flow theory (finite element simulation) and deformation theory 
is excellent. It is noted also that the singularity of the function cpi does not invalidate the results of 
the deformation theory model, as the deformation can be continued beyond the discontinuity line. 

In the second example, the following monotonic loading has been applied in two steps (see 
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Shear strain £12 Shear strain £12 

Figure 6: Response of the material in terms of stress-strain curves along the isochoric path AB of 
Fig. [5^: (a) shear stress (T12; (b) normal stress cJn. 


Fig. [7|: first a isochoric deformation (shearing in the plane X 1 -X 2 , from point O to point A in 
Fig. [7^), followed by a dilatational deformation (from point A to point B in Fig. [7^). This specific 
loading path has been chosen as it crosses the discontinuity line for the function (j )2 ■ The value of the 
function (1)2 along the dilatational path AB is shown in Fig. [7h: results of finite element simulation 
and deformation theory are in agreement. 




Figure 7: (a) Function (/)2 in the plane: elastic domain 02 = 0 (light grey), plastic domain 

with a discontinuity line separating the regions 02 > 0 (light blue) and 02 = 0 (light grey). The 
red dashed line indicates the loading path, (b) Function 02 along the dilatational path AB. The 
discontinuity is clearly visible and results of finite element simulation and deformation theory are 
in agreement. 


In Fig. IHl we present the stress-strain curves along the dilatational path AB of Fig. [7^. In particular, 
Fig. |8^ shows the shear stress <712 and Fig. [H]^ the normal stress an, both as a function of the 
dilatational strain en. Again, the agreement between flow theory (finite element simulation) and 
deformation theory is excellent. However, now the deformation cannot be continued beyond the 
discontinuity line, as the function 02 would become negative, which is unphysical. This is due to 
the fact that the Drucker-Prager model with associate plastic flow always predicts plastic dilatancy. 
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Normal strain £ii Normal strain en 


Figure 8 : Response of the material in terms of stress-strain curves along the dilatational path AB 
of Fig. [7^: (a) shear stress CJ 12 ; (b) normal stress an. 


3.4 Plane strain formulation of transmission conditions 

The two-dimensional plane strain formulation is relevant for the analysis if the conditions of the 
plane strain are satisfied in the system, i.e. the following components of the strain and stress tensors 
vanish within the interphase: 


£^■3 = 0, <7*3 = 0, j = 1,2,3, z = l,2. 


( 66 ) 


while the remaining components are functions of the independent variables x\ and X 2 only. 
Finally, the remaining component of the stress tensor can be evaluated from 


<733 — (^11 + £ 22 ) A (</> 1 , <;^ 2 )- 


(67) 


Taking into account 
invariants 


161 : 


and the Hooke’s law ([T]), we obtain the leading terms of the strain tensor 
du 2 CF22 


dx2 A-k2/i’ 


( 68 ) 


re = 1 + -[ — V = ^ ^22 , 1 ^12 

^ 3 Vdx 2 / 4 Vdx 2 / 3 (A+ 2/1)2 4 / 22 ' 


(69) 


Moreover, all the strain and stress components here are functions only of the spatial variable Xi. 
Finally, relationship (I16p remains valid, while (I17p is rewritten in the form: 


= 




12/l2 


16 /l 2 ' 


(70) 


Note that the transmission conditions dZD and dH]) are defined now for the two-components vectors 
^2 = [o' 2 i, o' 22 ]^, u = [ui, M 2 ]^ and the matrix A is given by 


r, _ 0 A 

\0 A -|- 2/1 j 


(71) 
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It should be noted that the transmission conditions for plane strain conditions can be written in 
the compact form as 

Icruj = 0, IfT22] = 0, ai2 = -Flduil, |u2l), 0-22 = F2{luil, |U21), (72) 

where the functions Fi and F 2 are given by 

Fi = F 2 =-^{X + 2jl){cj)i,4>2)lu2}. (73) 

4 Validation of the nonlinear transmission conditions for elasto- 
plastic pressure-dependent interphases 

The geometry of the considered structure is shown in Fig. |9l In the analysis the two bonded 
materials are assumed to be identical and linear elastic, with Young’s modulus E± = 72700 MPa 
and Poisson’s ratio i'± = 0.34. The geometrical dimensions are L = 10 mm, H = 1 mm, and 
2h = 0.01 mm. As a result, the small parameter takes the value e = 2h/H = 0.01, which is small 
enough to justify our asymptotic model. 




X+, Ai+ 

X2 : 



-C A 
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1 
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Figure 9: Geometry of the three-phase structure. 

A linear elastic-linear hardening Drucker-Prager model (pressure-dependent yielding) is assumed 
for the adhesive material within the interphase. The material parameters are defined in Tab. [T]and 
the corresponding stress-strain curve in uniaxial stress test is shown in Fig. |3^. For comparison 
with known results in the literature, a linear elastic-linear hardening von Mises model (pressure- 
independent yielding) is also considered. All the analyses performed with the von Mises model are 
in agreement with the literature and are not shown here. 

The commercial finite element code ABAQUS is used for the simulation of the mechanical behaviour 
of the three-phase structure and validation of transmission conditions (1721) . The transmission con¬ 
ditions are based on the deformation theory of plasticity developed in Section |31 whereas the FEM 
simulations are based on the more general theory of plastic flow. For this reason, only monotonic 
external loading is considered in the simulations. 

The two-dimensional FE-mesh is shown in Fig. [101 It has been generated automatically using the 
mesh generator available within ABAQUS with 4-node bilinear, reduced integration elements with 
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Material 

E 

(MPa) 

V 

Os 

(MPa) 

EP 

(MPa) 

uP 

a 

Drucker- 

Prager 

813 

0.3 

50 

81.3 

0.22 

0.1504 

von Mises 

813 

0.3 

50 

81.3 

0.53 

0 


Table 1: Material parameters used for the validation of transmission conditions. 


hourglass control (CPE4R). A strong mesh refinement is used within the intermediate elastoplastic 
layer. The mesh consists of 757774 nodes and 757599 elements. 

The validation of transmission conditions is performed as follows. From the FEM simulation 
we obtain the displacements at the top {x 2 = h) and bottom {x 2 = —h) of the intermediate layer 
and evaluate the jumps [ui] and [n2]- Then, we can compute in sequence: the strain invariants Jf 
and J| from (1161) and (1701) . the functions (pi and (p 2 involved in the deformation theory from (1651) 
and (I64p . the generalized Lame constants A and /2 from (125 p . and finally the tractions predicted by 
the transmission conditions Ti([[mi], [[n2]l) and T2 ([[mi]], [^ 2]) from (|73p . The final check consists in 
comparing the predicted tractions with those obtained from the FEM simulation. 



Figure 10: Two-dimensional FE-mesh. 


4.1 Tensile loading with fixed grips 

A monotonic tensile loading with fixed grips (ni(xi, 77/2) = 0, U 2 {xi, H/2) = V 2 ) is applied to the 
top of the bimaterial structure in the range of V 2 /H from 0% to 0.3% in 100 fixed increments, see 
Fig. [TT] 

The equivalent stress and the equivalent plastic strain along a horizontal line in the middle of 
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Figure 11: Tensile loading with fixed grips. 


the interphase are shown in Fig. [TJl every ten increments. It is noted that at the 10*^ {V 2 /H = 
0.03%), 20*^ {V 2 /H = 0.06%) and 30*^ {V 2 /H = 0.09%) increments the interphase is still in the 
elastic regime. A visible plastic zone appears in the middle of the interphase at the 40*^ increment 
{V 2 /H = 0.12%). At the 50*^ increment {V 2 /H = 0.15%), the interphase is fully plasticised. 




rci-coordinate [mm] 


xi-coordinate [mm] 


Figure 12: Distribution of equivalent stress (a) and equivalent plastic strain (b) along the interphase 
for different levels of deformation (tensile loading with fixed grips). 

Due to the symmetry of the considered problem (geometry and applied loading), three of the 
transmission conditions, i.e. [[CT12I = |o'22]l = 0 and Ti(|ui]], |tt2l) = < 712 , are identically satisfied, 
because |tti] = 0 and cri2 = 0 hold in this case. The remaining condition ^2(0, |u2l) = <722 has to 
be verified. 

In Fig. fTOl al. comparisons of the left and right hand sides of the condition ^2(0, |u2]l) = <722 are 
presented. The accuracy of the transmission condition is excellent, both in the elastic and plastic 
regimes, along the whole interface except for a very small edge zone, as shown in Fig. [T3l bl. where a 
magnification of the same functions in the interval 4.95 < xi < 5 at the right boundary is presented. 
The size of this edge region is typically 3-4 times the thickness of the interphase. 
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Figure 13: (a) Validation of the transmission condition ^2(0) [[^^2l) = o'22 for a pressure-sensitivity 
elastoplastic interphase (tensile loading with fixed grips), (b) Magnification of the same functions 
in the interval 4.95 < xi < 5 at the right boundary. 


4.2 Shear loading 


A monotonic shear loading /2) = m, U 2 {xi, H/2) = 

material structure in the range of Vi/H from 0% to 0.7% in 


0 ) is applied to the top of the bi- 
100 fixed increments, see Fig. [T4] 



Figure 14: Shear loading. 

The equivalent stress and the equivalent plastic strain along a horizontal line in the middle of 
the interphase are shown in Fig. [151 every ten increments. It is noted that at the 10*^ {vi/H = 
0.07%), 20*^^ {vi/H = 0.14%), and 30*^ {vi/H = 0.21%) increments the interphase is still in the 
elastic regime. A visible plastic zone appears in the middle of the interphase at the 40'’^ increment 
{vi/H = 0.28%). At the 50*'^ increment {vi/H = 0.35%), the interphase is fully plasticised. 

Due to the anti-symmetry of the considered problem (geometry and applied loading), two of the 
transmission conditions, i.e. [cri2] = [o'22]] = 0, are identically satisfied. The remaining conditions 
1^21) = cri2 and F2([ui]], |u2]) = (T22 have to be verified. In fact, due to the dilatancy 
of plastic deformation typical of the Drucker-Prager model, we have |tt2] 7^ 0 and, consequently, 
tensile stress develops within the interphase, <722 7^ 0. 

In Fig. [Tol a), comparisons of the left and right hand sides of the condition Ti(|ui], |u2l) = 0'i2 
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Figure 15: Distribution of equivalent stress (a) and equivalent plastic strain (b) along the interphase 
for different levels of deformation (shear loading). 


are presented. The accuracy of the transmission condition is excellent, both in the elastic and 
plastic regimes, along the whole interface except for a very small edge zone, as shown in Fig. HHKb), 
where a magnification of the same functions in the interval 4.95 < xi < 5 at the right boundary is 
presented. The size of this edge region is typically 3-4 times the thickness of the interphase. 



Figure 16: (a) Validation of the transmission condition Ti([[tii], [[n2l) = 0'12 for a pressure-sensitivity 
elastoplastic interphase (shear loading), (b) Magnification of the same functions in the interval 
4.95 < xi < 5 at the right boundary. 

In Fig. fTW a). comparisons of the left and right hand sides of the condition F2(|ui], |u2]l) = (T22 
are presented. In this case, the accuracy of the transmission condition is still excellent in the elastic 
regime {vi/H = 0.14%) and in the early stage of plastic deformation {vi/H = 0.28%). However, 
at later stages of plastic deformation {vi/H = 0.7%) the prediction of the transmission condition 
appears to be less accurate. 

We observe that, the ellipticity conditions for the generalized elastic constants are not satisfied 
in the considered stages of plastic deformation {vi/H = 0.28% and V\/H = 0.7%), see Fig. [181 In 
particular, the generalized bulk modulus K becomes negative and, consequently, the generalized 
Poisson’s coefficient exceeds 0.5. This behaviour is intrinsic of the Drucker-Prager model with 
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Figure 17: (a) Validation of the transmission condition F2([[^*i]) = <722 for a pressure-sensitivity 

elastoplastic interphase (shear loading), (b) Magnification of the same functions in the interval 
4.5 < xi < 5 at the right boundary. 


associated plastic flow, which always predicts plastic dilatancy. 




xi-coordinate [mm] 



Figure 18: Generalized elastic constants fl {a), K (b), u (c), for different levels of deformation (shear 
loading). 


4.3 Combined loading 

A combined loading is applied to the top of the specimen in such a way that the same displacements 
are prescribed in the xi as well as the X 2 directions {ui{xi,H/2) = U 2 {xi, H/2) = v) in the range 
of v/H from 0% to 0.35% in 100 fixed increments, see Fig. [191 

The equivalent stress and the equivalent plastic strain along a horizontal line in the middle of the 
interphase are shown in Fig. [20l every ten increments. It is noted that at the 10*^ {vi/H = 0.06%) 
and 20*'^ {vi/H = 0.12%) increments the interphase is still in the elastic regime. A visible plastic 
zone appears in the middle of the interphase at the 30’’^ increment {vi/H = 0.18%). At the 40**^ 
increment {vi/H = 0.24%), the interphase is fully plasticised. 

For this type of loading, the two conditions Fi([[ui]], |u2]l) = <7i2 and F2([[ui]], |tt2]]) = <722 have to be 
verified. In Fig. l2lT ai. comparisons of the left and right hand sides of the condition F2([ui], [[n2]l) = 
(722 are presented. The accuracy of the transmission condition is excellent for all stages of deforma¬ 
tion. In Fig. l2^ al. comparisons of the left and right hand sides of the condition Ti(|ui], [[U2I) = <7i2 
are presented. 


20 













































1 t t t t ' 

f t t f f 1 



A 

-► Xi 

1 

1 


77777Z77777Z77777Z77777777777777777777777777777777777 


Figure 19: Combined loading. 
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Figure 20: Distribution of equivalent stress (a) and equivalent plastic strain (b) along the interphase 
for different levels of deformation (combined loading). 



Figure 21: (a) Validation of the transmission condition F2([[^*i]) = <722 for a pressure-sensitivity 

elastoplastic interphase (combined loading), (b) Magnification of the same functions in the interval 
4.95 < xi < 5 at the right boundary. 
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Figure 22: (a) Validation of the transmission condition Fi([[rti], [[U2I) = <712 for a pressure-sensitivity 
elastoplastic interphase (combined loading).(b) Magnification of the same functions in the interval 
4.95 < xi < 5 at the right boundary. 


5 Conclusions 


We have proven that a thin interphase consisting of a soft elasto-plastic pressure dependent material 
can be effectively modelled by an imperfect nonlinear interface. The latter is governed by nonlinear 
transmission conditions prescribed along the infinitesimal interface separating two different materials 
and incorporating the material properties of the adhesive material. The transmission conditions 
evaluated here differ essentially from those previously delivered for the case of von Moses plastic law 
13|, ll4|- Indeed, the phenomenological functions (f)i and ^2 describing the equivalent elasto-plastic 


interface are not smooth and exhibit a singular behaviour having clear physical sense. Moreover, 
the conditions evaluated here are universal, in comparison with those delivered previously |l6l |. 
and can be used for arbitrary stress-state of the considered three-phase body. This makes the 
evaluated transmission conditions widely applicable and ready to use in FEM analyses (for instance 
as traction-separation laws in cohesive element formulations) in order to model very thin adhesive 
layers. This approach avoids extreme mesh refinements usually needed in standard approaches. 

The transmission conditions are validated by numerical examples based on accurate finite el¬ 
ement simulations. The accuracy and efficiency of the proposed approach prove to be very high 
for different monotonic loading conditions. However, as it is the case of any imperfect interface, 
the conditions are valid along the whole interface apart from a small region close to the sample 
boundaries where the edge effect is playing a decisive role. To tackle this phenomenon, one needs 
additionally to use special FEM elements properly modelling the local stress-strain field. Such 
strategy (imperfect interface approach enriched by the edge special elements) is an efficient way 
to compute the details of a structure containing extremely thin soft elasto-plastic interphases of 
various plastic properties, at least in case of monotonic loading. In case of cyclic load, the range of 
applicability of the transmission conditions is still to be analysed. 
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